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Recently a series of indications have been put forward suggesting the presence of two 7-ray lines 
at 110-130 GeV (centered at 111 and 129 GeV). Signals of these lines have been observed toward 
the Galactic center, at some galaxy clusters and among some of the unassociated point sources of 
the 2 years Fermi catalogue. Such a combination of signals could be generated by dark matter 
annihilations in the main dark matter halo, its substructures and nearby galaxy clusters. We 
discuss here the consistency between the number of events observed at the line energies in the sky 
and the predictions using results from the Via Lactea II numerical simulation and extrapolations 
below its mass resolution, taking into account that the annihilation cross-section to the lines can be 
estimated from the Galactic center signal. We find that some extrapolations to small substructures 
can naturally account for the point sources signal, although the hypothesis of background only cannot 
be rejected. We also study the morphology of the 7-ray sky at the 2 lines energies, testing different 
Galactic diffuse background models to account for interstellar medium uncertainties, and different 
assumptions on the DM diffuse component profile. We find from template fits that within reasonable 
diffuse background uncertainties the presence of a spherical halo component is preferred and that 
cuspier dark matter halo profiles are also preferable even from the full sky fit. We finally check the 
impact of a dark disk component, suggested by cosmological simulations that include baryons and 
find that thin dark disks can not be disfavored, thus possibly accounting for the preferentially closer 
to the Galactic disk distribution of the point sources lines signal. 
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I. INTRODUCTION 

The possible identification by the Fermi Gamma- 
ray Space Telescope of a signal compatible with the 
monochromatic photon emission due to pair annihilations 
of cold dark matter (DM) particles has recently been one 
of the most debated topics. Originally, [l[ has suggested 
the detection of a line at 129.8 ± 2A~^[ 3 GeV with a 3.3a 
significance in a wide window toward the Galactic cen- 
ter (GC). A similar signal has been indicated by @ at 
127.0 ± 2.0 GeV with 5. Oct significance. A pair of lines 
with energies of 110.8 ± 4.4 and 128.8 ± 2.7 GeV can 
alternatively explain 7-ray excess with 5.4cr significance 
0. Similarly the line signal at ~ 130 GeV has also been 
found by [j| with suggesting the presence of 2 lines 
at ~ 110 and 130 GeV. Both line signals are in agree- 
ment with constraints from line searches of the Fermi 
collaboration Q and indicate a preference for dark mat- 
ter (DM) annihilation rather than decay 0, H, H, H • The 
Fermi collaboration, motivated by the results of [H-Q, 
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has performed an analysis oriented toward the GC and 
some of those results are presented in [13, [ll[ . Similar 
to [H, [HI , no obvious systematic error has been found 
to account for the amplitude of the line signal measured 
by [E (see though [13] )■ The Fermi results confirm 
a line-like signal at E~130 GeV at 4cr, or E~135 GeV 
at 3.3c after reprocessing the data to take into account 
the shift of the reconstructed energy with time. Yet, 
some part of the amplitude may be related to limb pho- 
tons [lOj . Adding information on the performance of the 
instrument's energy reconstruction decreases the signif- 
icance of the signal Thus a conclusive answer on 
whether the line signal is a systematic error identified 
toward the GC or a signal of DM annihilations, has not 
been provided yet. 

A monochromatic gamma-ray flux is expected in most 
scenarios in which DM is in the form of weakly interacting 
massive particles (WIMPs), since two-body annihilation 
final states containing photons arise at the 1-loop level. 
At the same time, in such a framework, it is foreseen 
that tree-level WIMP annihilations into other SM final 
states, in turn hadronizing and/or decaying into p, p , 
us and 7s, would give sizable yields with continuum en- 
ergy spectrum on top of the monochromatic 7-ray yield. 
Yet, no clear indication of 7-ray excess, other than the 
lines, has been found toward the inner few degrees of 
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the GC, leading to the extraction of constraints on the 
continuum [El. ttWtYi\ . 

Another aspect of a WIMP annihilation signal is that 
one should expect to see line signals at the same energy 
and annihilation cross-section toward other dark mat- 
ter targets. In [18], a 130 GeV line signal toward known 
galaxy clusters has been suggested. There have been also 
indications for two lines at 111 and 129 GeV in unasso- 
ciated point sources which would imply the DM annihi- 
lation in substructures 0, H3] [§l[ (see also [H[ for an 
alternative interpretation) . 

We assume that the line signal from the unassociated 
point sources is indeed of DM origin with the same anni- 
hilation cross-section to the lines as is estimated from the 
GC. We then confront that signal with predictions from 
cosmological simulations such as Via Lactea [22-24J . In 
section |TT1 we discuss the 7-ray data that we use and our 
general assumptions for the background and the DM den- 
sity distribution in the Galaxy. As a general reference we 
use the total number of 7-ray events with energy 111 ± 5 
GeV and 129 ± 6 GeV. These events are taken to be of 
both DM and diffuse/point source 7-ray background ori- 
gin. The comparison of the observed 7-ray data with 
the predicted contribution from substructures within the 
Galaxy's virial radius is done in section IIIII Our aim is 
to conclude on whether the line signal of point sources 
can be physically associated to the same energy line sig- 
nal toward the GC and under what assumptions on the 
substructure distribution. On the CDM simulation side 
we use VLII subhalo distribution data [2J] and also ex- 
trapolate the VLII simulation mass function to smaller 
subhalo masses. 

The observed non-isotropic distribution of DM line(s) 
point sources in the sky could be explained by the pres- 
ence of a strong dark disk. Such a dark disk would also 
have an impact on the diffuse distribution of the DM orig- 
inated line photons. To study the diffuse 7-ray sky at the 
energies of the 2 lines, one needs predictions on both the 
Galactic diffuse backgrounds and on the DM diffuse con- 
tribution. In section ITV1 we test that possible contribu- 
tion to the two 7-ray lines from DM annihilations in the 
Galaxy and diffuse backgrounds by doing a template fit. 
The importance of template fits is that one can take into 
account the different morphologies of the various diffuse 
components. We study the impact on the significance of 
a DM signal on the Att sky (including the GC) of dif- 
ferent assumptions for the Galactic diffuse background, 
related to physical properties of the interstellar medium. 
We also test different assumptions (and thus different 
templates) on the main spherical halo density profile, on 
the significance of the dark disk component to the local 
dark matter density and its thickness and finally on the 
contribution of the dimmer DM subhalos that would also 
add to the diffuse 7-ray sky. We consider that the main 
contribution from the brighter DM structures have been 
already observed by [19J and exclude them from the 7-ray 
fits. We also derive upper limits on the diffuse emission 
from annihilations in the main DM halo, and give our 



conclusions in section IVl 



II. GAMMA-RAY DATA, DIFFUSE 
BACKGROUND AND DM DISTRIBUTION 
ASSUMPTIONS 

The Fermi Large Area Telescope publicly available 
events are categorized in different classes based on the 
expected level of cosmic ray (CR) contamination. In this 
work we use the ULTRACLEAN events class which is 
the cleanest 7-ray events sample. There are 686(744) 
photons with energy between 111 ± 5 GeV and 611(668) 
photons with energy between 129 ± 6 GeV in the 4 yr 
(4.4 yr) full sky Fermi ULTRACLEAN class 7-ray data, 
with the quoted energy ranges representing the relevant 
energy dispersion for these lines 1621 . In section [TTTI where 
we compare with the findings of [19( , we use the 4yr sam- 
ple since it approximates their events sample, while in 
section QV] we use the slightly longer period of 4.4 yr . 

The emission of the diffuse Galactic 7-ray background 
above 100 GeV is dominated by the ir° contribution, i.e. 
the decay of mesons produced by inelastic collisions of 
CRs with the interstellar medium (ISM) gas, and by 
the up-scattering of low energy photons of the interstel- 
lar radiation field (ISRF) from high energy CR electrons 
(inverse Compton scattering). The morphology of these 
components on the sky is different mainly because of dif- 
ferent morphologies of the ISM gas density and the ISRF 
energy density in the Galaxy. Moreover the energy loss 
of CR electrons and protons during their propagation in 
the Galaxy is different. The bremsstrahlung radiation of 
CR electrons at these energies is completely subdominant 
but it is included in our code. 

To compute the diffuse 7-ray background, we use the 
DRAGON package [H^l with a new ISM gas model [Hi 
that ensures good agreement with 7-ray spectral data be- 
tween 1 and 200 GeV in the full sky and subsections of it 
[29l ]. We ignore the contribution of the "dark gas" (not 
related to DM substructures) whose uncertainties are sig- 
nificant in the inner 5° in latitude [3(| [H|. Based on the 
relevant uncertainties (see [29]), we allow for different 
assumptions on the ISM gas and the ISRF which influ- 
ence the 7T° and the inverse Compton 7-ray emissivities 
respectively. 

In the case where there are 2 lines as has been indicated 
by @,|1], the energy of these lines is centered at 128.8±2.7 
and 110.8 ±4.4. GeV 0. The lines come from either the 
combination of 27&Z7 lines or from the Z7&/17 lines. 

In [17j , five individual modes/channels of DM anni- 
hilation: XX — * W + W~ , XX — ► bb, XX — ► t + t~, 
XX — ► an d XX — > e + e~ have been studied. The 

limits on the DM annihilation cross-sections based on 
their contribution to the continuum 7-rays spectrum in 
the I Z |< 5°, | 6 |< 5° observation window have been de- 
rived. Typically, DM models have sizable branching ra- 
tios into more than one of these channels. Yet apart from 
the XX — ^ channel and mainly the XX — ^ e + e~ 
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channel, in all the other annihilation channels to SM par- 
ticles with a continuum spectrum, the 7-ray DM signal 
at 111 and 129 GeV can not be explained/mimicked by 
the continuum spectrum. Thus it originates from the 
annihilation into Z7 and 27. For XX — * e + e~ and 
XX — > P + P ~ 1 the final state radiation (FSR) and virtual 
internal bremsstrahlungJVIB) can contribute to the line 
signal as discussed in pjl HH . 

For simplicity we assume that the DM induced 7-rays 
with energy 111 ± 5 GeV and 129 ± 6 GeV come from 
the annihilation of a 129 GeV DM particle into Zj and 
27 respectively. Alternatively, these 7-ray lines could 
come from /17 and Zj for the case of a 142 GeV DM 
particle. The relevant ratio of the luminosity of two lines 
is taken to be 1/2 for the 111/129 GeV lines as suggested 
in (l9fl . thus, for the case of 129 GeV DM particle, the 
annihilation cross-sections to Zj and 27 are assumed to 
be the same. 

For the DM distribution we assume that it is a combi- 
nation of a spherically symmetric " main" DM halo and a 
dark disk (DD). For the main halo we assume a spherical 
Einasto DM profile: 



Ps P h(r) = PEm exp < -- 



-1-1 

r, 



(1) 



using S = 0.13,0.17,0.22 [33j with r c = 20 kpc. The 
values of 5 = 0.13(0.22) result in a more (less) cuspy DM 
distribution. The density normalization parameter pEin 
is set in terms of the local DM density, after including a 
contribution of the DD. 

The profile of the DD component is assumed to be 
described by [34] 



Pdd{R,z) = po DD exp 



1.68 (R Q - R) 



Ri 



/2 



exp 



0.693 \z\ 

Zl/2 



(2) 

where R1/2 and Z1/2 are the half mass scale lengths in the 
Galactic plane and perpendicular to the Galactic plane, 
respectively, and Rq = 8.5 kpc. Here R is the cylindrical 
radial coordinate. 

The ratio of the local DM density in the dark disk to 
the local DM density in the spherical halo po DD / Po aph 
typically range between 0.2-1.5 [341 ], with the higher ra- 
tios being related to higher mass densities in the thick 
stellar disk rather than in the thin stellar disk. The 
thick stellar disk can be populated by thin stellar disk 
stars, if the thin stellar disk gets heated by very massive, 
high-redshift mergers. Another cause could be multiple 
pro-grate and low inclination mergers [34| . 

In the template analysis performed below we will re- 
strict to the case: 

a/2 = po DD /(po sph +Po D d) < 0-5. (3) 

fixing [H[Ii: 

Po sph + POnn = 0.4 GeV cm" 3 . (4) 



Regarding the dark disc thickness, some authors [37| 
have suggested thicker disks, while thinner and less sig- 
nificant dark disks can also be the case; keeping in the 
parametrization of eq. H]i?i/2 = 11-7 fixed 34j . we will 
test the half mass scale length values of zx/i = 0.5, 1.0, 
1.5, 3.0 kpc. 

In the standard model for cosmology, cold DM struc- 
tures form hierarchically, with small DM halos collaps- 
ing first and subsequently merging into larger and larger 
objects. Since tidal disruption may only be partially ef- 
fective, massive DM halos, such as the halo of our own 
Galaxy, are expected to contain a vast population of sub- 
halos, with mass spanning from a tiny seed mass up to a 
fraction of the hosting halo mass. The minimum mass is 
essentially associated to the free-streaming scale of DM 
particles, in turns depending on their temperature of ki- 
netic decoupling in the early Universe; for WIMPs the 
minimum mass can be as small as about m cu t = 10 _6 Mq 



[38l l39jj , much lighter than the dwarf galaxy scale, pos- 
sibly to the smallest environment which can host stellar 
populations and hence a luminous counterpart. Because 
of the highly non-linear nature of the merging process, 
up to now the only efficient technique to model in detail 
DM halos is the use of numerical N-body simulations; in- 
deed large populations of substructures have been found 
in such studies. We will assume as primary reference in 
our analysis results from Via Lactea II (VLII)[4(|, one 
of the highest resolution simulations up to date of Milky 
Way-sized CDM halo (virial mass M h = 1.9 x 10 12 M Q ), 
with over one billion DM "particles" and nominal mass 
resolution of about 4100 M Q (numerical effects appear 
to enter well above this scale, possibly affecting the sub- 
halo mass spectrum up to about ~ 3 x 10 6 Af Q ). In 
our analysis we will discuss both the DM pair annihi- 
lation associated to individual DM substructures as well 
as the collective effect from the whole subhalo popula- 
tion; in both respects, the resolution of the simulations 
appears insufficient to properly model the expected sig- 
nals. Our approach will then be to use the simulation 
results to properly calibrate the necessary extrapolations 
to smaller masses: tuning, at given Galactocentric ra- 
dius, the subhalo pericenter distribution and applying a 
recipe for taking into account tidal stripping effects, we 
derive a model which reproduces fairly well the subhalo 
mass function and the distribution in halo concentration 
as a function of radius in the VLII simulation, and we 
use it as a prediction below its resolution (some details 
about our approach are given in Appendix [AJ . 

The general trends in the DM subhalo distribution can 
be understood from the fact that more massive objects 
are more prone to tidal stripping than the less massive 
ones, because they typically have smaller average density, 
reflecting the fact that they collapsed later in the cosmic 
history at a lower average background density. As the 
result, when going toward the center of the host halo, 
the average subhalo density increases and the average 
mass decreases. Also, in the inner part of a DM halo the 
tidal forces become stronger, possibly making the sub- 
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FIG. 1: Average DM squared density (p 2 ) associated to sub- 
halo population, as a function of Galactocentric radius r g , 
for the case of a biased distribution (solid lines) and an anti- 
biased distribution (dashed lines) , and for a few assumptions 
on the values for subhalo mass function a and minimum sub- 
halo mass m cu t (see the text for details); also shown is the 
density squared of the smooth halo component. 

halos spatial distribution to be anti-biased with respect 
to host density profile: this tendency has been found in 
some numerical simulations, see, e.g., 0, while it 
has also been claimed that biased or anti-biased distri- 
butions may just stem from selection effects [23], |4(j . The 
model in our extrapolation from the VLII results gives a 
biased distribution; to bracket uncertainties we consider 
also an anti-biased distribution which is instead taken 
from [241 ] . The other parameters entering most critically 
in our analysis are the spectral index a for the subhalo 
mass function and the minimum subhalo mass m cut . In 
Fig. [T] we plot, as a function of Galactocentric radius, 
the average DM squared density associated to the full 
subhalo population, for two choices of the spectral index, 
i.e. a = 1.9 and a = 2, three sample values of the mini- 
mum subhalo mass, and for the biased (solid lines) and 
anti-biased (dashed lines) distributions; as a comparison 
the density squared of the smooth DM halo component, 
Eq. [TJ with S — 0.17, is also given. 



III. DM SUBSTRUCTURES IN THE MILKY 
WAY AND THE 2 7-RAY LINES 

In [l||, 16 unassociated point sources have been identi- 
fied with at least one 100-140 GeV photon for SOURCE 
event class . Among those photons, there are 9 photons 
in the ULTRACLEAN event class sample with energies 
~111 and ~129 GeV. Those photons from unassociated 
point sources may imply a signal of DM annihilation in 
close by substructures. That ansatz can be compared 
with predictions from cosmological simulations. 

The 2 yr point source catalogue [42j has 575 unas- 
sociated point sources. In [HJ|, the 9 ULTRACLEAN 



photons at ~ 111 and 129 GeV lie within a 0.15°/0.3° 
radius for FRONT/BACK converted events covering an 
area of 0.07/0.28 square degrees. Thus the 575 unasso- 
ciated point sources cover at most 1.6 x 10 2 deg2 out of 
4.1 x 10 4 deg2 and would give a conservative upper es- 
timate of 5.06 photons out of the 1.3 x 10 3 63]. The 
probability that 9 or more ULTRACLEAN class photons 
out of the 1.3 x 10 3 photons of the 4-7T sky fall within 
the area covered by the 575 unassociated point sources is 
p = 0.0721 I64J. Yet there may also be a bias toward dis- 
covering point sources around single high energy events 
[43| . While by itself not a strong statistical deviation, 
the coincidence in energy with the GC signal inclines us 
to test for alternative possibilities and whether some of 
these photons come from unassociated point sources that 
are DM subhalos. 

The number of photons that we receive from a single 
subhalo with luminosity L = J p 2 sub dV and line of sight 
(los) distance A from us, for channel ch — 77 or 7Z is 
given by: 

N ch _ Kfch ( a Av)ch L T exp A exp 

7 2 m 2 4ttA 2 ' { 1 

where m x = 129 GeV is the DM particle mass, and 
N^ h — 1(2) for 7^(77). r exp and A exp are the detec- 
tor's exposure time and effective area for photon's en- 
ergy of 129 GeV, respectively. In this work, we use 
for the averaged Fermi-hAT exposure after ~ 4 years 
r exp x A exp = 1.22 x 10 11 cm 2 s. For this section, we con- 
sider two values of annihilation rates: (o"i>) 77 = (av)z-y — 
0.98 x 10~ 27 cm 3 s" 1 , a value derived assuming our de- 
fault smooth component DM density profile (no DD) and 
fitting the monochromatic signal in the region | I \ < 5° & 

1 b |< 5° [13; and {ov). n = {av) Zl = 3 x 10" 28 cm 3 s"\ 
which fits better the whole sky region (see section IIVI) . 

To quantify the possible impact of substructures in the 
line photons on the sky we ask the following questions: 

A) How many subhalos give 2 or more photons in to the 

2 7-ray lines? 

B) How many photons (in the two lines energies) do we 
get from all the subhalos that give a more than 0.1 pho- 
tons? 

The difference between the number of photons from the 
entire subhalo population received and the answer to 
question (B) is a proxy for the diffuse gamma-ray flux 
to the two lines from DM substructures gravitationally 
bound in the main DM halo. We will refer to these pho- 
tons as " DM substructure diffuse" . 

In Table Q] we consider first the single subhalo sample 
from the VLII simulation [40( (no extrapolation below the 
mass resolution at this level) and compute answers to the 
questions formulated above, averaging over results ob- 
tained for 100 random choices for the position of the ob- 
server, all at fixed Galactocentric distance Rq — 8.5 kpc 
(the average is performed to wipe out fluctuations in- 
volving effects of nearby subhalos or voids in a single 
random choice). We then turn our analysis to the ex- 
trapolated subhalo populations focussing on the biased 
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VLII 


(0) 


0.213 (0.024) 


biased - case I 


0.0198 (0.00344) 


0.473 (0.0874) 


biased - case II 


0.0139 (0.0024) 


0.342 (0.0618) 


anti-biased - case I 


0.0746 (0.0176) 


1.24 (0.296) 


anti-biased - case II 


0.0898 (0.0196) 


1.62 (0.361) 



TABLE I: Relevance of substructures for detection of the 
monochromatic photons as referred to the questions A to B 
as posed in the text, using cross section which fits GC (fits 
the whole sky). Answers are provided in case of the subhalo 
sample from the VLII simulation itself and from extrapola- 
tions of it in case of biased or anti-biased distributions with 
the parameter choice m cu t = 10 3 Mq & a = 1.9 (case I) and 
iTi cu t = 10 -6 Mq & a = 2 (case II). Changing the overall 
normalization of the subhalo number density would shift the 
results provided in the table accordingly, e.g., by a factor of 
about 2 if adopting the normalization of the Aquarius simu- 
lation [m. 



and anti-biased distribution and considering two possible 
extrapolations for lowest subhalo mass and subhalo mass 
spectral index, bracketing extreme possibilities, case I 
corresponding to (m cutl a) = (10 3 M©, 1.9) and case II to 
(m cut ,a) = (1CT 6 M Q , 2). The number of line photons 
received is computed including only subhalos within the 
virial radius since we find that for VLII subhalos only 
0.5% of the total photon comes from substructures lying 
outside it. On the other hand there is a further uncer- 
tainty one should be careful about: by tuning our subhalo 
model to the VLII results we are fixing the normaliza- 
tion of subhalo number density (above mass resolution) 
for the Milky Way halo according to that specific real- 
ization; this is a quantity which actually has a certain 
scattering among different halo realizations and different 
simulations, e. g., referring to the results of the Aquar- 
ius simulation [4l| we should increase such normalizing 
of about a factor of 2, shifting results in Table fl] by the 
same factor. 

As a further test, in Table ILT1 we report the total num- 
ber of monochromatic photons expected from our entire 
subhalo populations under a set of different assumptions 
for spectral index and cutoff mass. The results are again 
shown for biased and anti-biased distributions. 

If the photons from unassociated point sources are 
from DM annihilation in substructures, their number (9 
ULTRACLEAN events) probes the number of photons 
from the brighter substructure subsample. Considering 
that the number of photons originating from subhalos 
that emit more than 0.1 photons in the two lines (i.e. 
the results of Question B in Table |T]), is indicative of 
such number of photons, we compare the probability of 
having observed 9 (or more) photons for a DM signal 
calculated in the VLII sample (0.213 photons), in ex- 
trapolation for biased distribution and case I (0.473 pho- 
tons), for anti-biased distribution and case I (1.24 pho- 
tons), for biased distribution and case II (0.342 photons), 
and for anti-biased distribution and case II (1.62 pho- 



tons). These probabilities are p = 0.0874 for the back- 
ground plus the DM signal stemming from VLII sam- 
ple, p = 0.108(0.185) for background plus the DM sig- 
nal in the extrapolation for biased (anti-biased) distri- 
bution in case I and p — 0.0975(0.23) for background 
plus the DM signal in the extrapolation for biased (anti- 
biased) distribution in case II. Using the normalization 
from Aquarius simulation would increase the number of 
subhalos and received photons by a factor of 2, shift- 
ing the probabilities to, respectively, p = 0.153(0.344) 
and p = 0.127(0.449). Thus the most conservative VLII 
assumptions case is marginally favorable than the just 
background case. When extrapolating below the mass 
resolution, probabilities increase further, reaching rele- 
vant levels in optimistic extrapolations. On the other 
hand, using the cross section which fits the whole sky 
(see section HV|) . we don't see much differences in the p- 
values from having just a background signal. Also, from 
Question A, there are no subhalos expected to give more 
than 2 photons. One must keep in mind that lowering 
the value of cross section by a factor of 3, by going from 
av which fits GC better to the one which fits the whole 
sky better, does not simply reduce all the values in Ta- 
ble I by the same factor. This is because the photons 
produced by each subhalo will decrease; so that some 
subhalos which previously gave more than, say, 0.1 pho- 
tons; will now give less. The photons coming from such 
subhalos are not included anymore. (Similarly with the 
numbers of subhalos which give more than 2 photons.) 
However, in Table II, all values do lower by the same 
factor, because they are the number of photons coming 
from all subhalos. 

The differences between the numbers of photons that 
originate from all DM subhalos (Table [TT]) and the num- 
bers of photons that originate from DM subhalos that 
contribute 10" 1 lines photons or more (Q.B) are conser- 
vative probes to the diffuse contribution from the DM 
subhalos at ~ 111 and 129 GeV. The VLII sample gives 
1.21 — 0.213 = 1.0 DM substructure diffuse component 
photon, extrapolation for the biased (anti-biased) distri- 
bution in the easel 5.46-0.473 = 4.99 (3.9-1.24 = 2.66) 
photons, while in case II 96 - 0.342 = 95.7 (87 - 1.62 = 
85.4) photons. An upper (rough) limit to the DM sub- 
structure diffuse ~ 111 and 129 GeV photon component 
can be derived by considering it approximately isotropic 
and then counting the ~ 111 and 129 GeV photons laying 
above b |> 60° times 7.46 (the ratio of 47r to the area of 
the sky with | b |> 60°). There are 40 111±5 GeV and 30 
129±6 GeV photons above | b |> 60°, i.e. an upper esti- 
mate of the ~ 111 and 129 GeV photons in the isotropic 
component is 522 photons; thus significantly larger than 
the 1.0, 4.99 (2.66), or 95.7 (85.4) DM substructure dif- 
fuse component photons predictions from VLII and bi- 
ased (anti-biased) case I and II. In section llVl we have a 
more model-dependent estimate of the isotropic 111 and 
129 GeV 7-rays component, which though decreases the 
isotropic component photons down to ~190-230. 

While for the anti-biased distribution we find more 
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Index "a" 




biased 


anti-biased 


2.0 


l.Ox 10 _b 


96 


87 


2.0 


1.0 


20.8 


20.4 


1.9 


1.0xlO~ b 


16.3 


10.2 


1.9 


l.OxlO 11 


5.46 


3.90 


1.9 


2.0xl0 4 


4.02 


2.99 



TABLE II: Number of 111 and 129 GeV lines photons con- 
tributing to the DM subhalo diffuse lines component, for vari- 
ous choices of subhalo distributions, using cross section which 
fits GC (using cross section which fits whole sky instead, will 
scale all values by a factor of ~ 3/9.8 = 0.3). We show re- 
sults for different subhalo mass function spectral index a and 
lower mass cut-off m cut , and for biased and anti-biased distri- 
butions. 

events from fewer sources than the biased distribution, 
since the subhalo concentration has a sharper dependence 
on Galactocentric radius with higher luminous subhalos 
closer to the GC, here the trend is reversed when sum- 
ming over the whole population of dim sources. 



IV. DIFFUSE 7- RAY LINES EMISSION FROM 
DARK MATTER ANNIHILATION 

As discussed in [IT] the 7-ray lines centered at 111 and 
129 GeV that are observed in the sky originate from com- 
bination of sources. We mask out the contribution from 
the 16 point sources detected at [l9( with one photon at 
energies between 100-140 GeV. We also mask out the ex- 
tended sources (galaxy clusters) where a similar excess 
of 100-140 GeV 7-rays has been observed [l8|. We use 
a mask of 0.5° in radius for each of the 16 the point 
sources of [19|] and a 4° radius mask for the targets of 
[iljl ]. The remaining contribution to the 111 and the 
129 GeV lines may come from the diffuse 7-rays emis- 
sion due to DM annihilations in the main halo and its 
dark disk, from background 7-rays produced in the Milky 
Way, from DM annihilation at small scale substructures 
in the Milky Way (that we have not yet identified as 
point sources) and from other isotropically distributed 
extragalactic astrophysical sources. DM annihilation in 
extragalactic structures and CR contamination will give 
an additional isotropic component. In Fig. [5] we show 
the 4-7T sky after implementing our mask on the 16 point 
sources and the 6 extended ones. 

The DM annihilation rate in any part of the Galaxy is 
given by: 

r Q „ ti = - — ?(<j ann I v |) (6) 

2mJ 

x (plph + PDD + 2 Psph ■ PDD + Plub) i 

with (a ann \v\) the annihilation cross-section taken to be 
the same for both the DM particles in the dark disk, 
the spherical halo and the substructures. For the case of 
Sommerfeld enhancement these cross-sections are in gen- 



Fermi-LAT 4.4yr 3 Masked Map 




FIG. 2: 7-ray events (ULTRACLEAN class) with energy of 
111±5 and 129±6 GeV after 4 yrs of collection by Fermi-LAT. 
We mask out the 16 point sources of [19J and the 6 extended 
sources of [3 (see text for more details). We present 7-ray 
events in Mollweide projection using HEALPix [44 ] . 

eral different, because of the dependence of the annihila- 
tion cross-section to the velocity dispersion of dark mat- 
ter (45rl48| and the fact that for the DM particles in the 
DD the dispersion is suppressed by a factor of 5-6 com- 
pared to that in the spherical halo components [34] ; with 
subhalos having even lower velocity dispersions. Thus for 
Sommerfeld enhanced models the dark disk contribution 
to CRs and 7-rays can be much more significant 49] (see 
also |50j for a discussion on the impact of subhalos). 

In fitting to the 2 7-ray line full sky data we probe the 
prompt 7-ray DM annihilation component of the spec- 
trum which is directly related to the annihilation rate in 
eq. [6] 

We use the masked full sky data with energies 1 1 1 ± 5 
GeV and 129±6 GeV. We perform a maximum likelihood 
fit calculating the log-likelihood based on [5lj |: 

In £ = 2J h In fii - fii - ln(fc l !) , (7) 

i 

where fii is the model of linear combination of templates 
at pixel i, and k is the map of observed counts which is 
just the single 111 ± 5 GeV and 129 ± 6 GeV 7-ray Fermi 
masked map. Our diffuse 7-ray model is composed of 6 
templates with 4 free parameters: 

Hi = N ■ Back; + A • [(2 - a) 2 ■ SphDM^ (8) 
+ a 2 ■ DarkDiski + a(2 - a) ■ MixedDM, 
+ SubDM,;] + B ■ ISO*. 

The Backi template comes from our DRAGON run and 
is kept fixed modulo a normalization N for a specific 
set of assumptions on the ISM and ISRF, the SphDM,; 
refers to the term in cq. [5] that is proportional to p 2 ph - 
The DarkDisk^ template refers to the term in eq. [5] pro- 
portional to p 2 DD and the MixedDM^ template to the 
2p S p/i • Pdd term. The contribution of dim Galactic DM 
subhalos to the diffuse 7-rays is included in the " SubDM" 
(related to the p 2 ub term in eq.O, where as described ear- 
lier we have masked out the brightest possibly detected 
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FIG. 3: The 7-ray templates as can be used in eq. [8] Top left: the Galactic diffuse background template including the 7r°, 
inverse Compton scattering and bremsstrahlung components at energies of 111±5 and 129±6 GeV. We show the "Back A" 
model. Top right: the diffuse DM spherical halo component assuming an Einasto with 5 = 0.13 profile. Bottom left: diffuse 
emission from the combined spherical and dark disk DM distributions, assuming a maximal dark disk (a = 1), which results in 
showing the "SphDM"+ "DarkDisk" + "MixedDM" combined template. We use Einasto profile with 8 = 0.13 for the spherical 
and Z1/2 = 0.5 kpc for the scale hight of the dark disk. Bottom right: the DM diffuse subhalo template "SubDM" for biased 
distribution with m cu t = 1O _6 M0 together with the spherical DM halo "SphDM" template (Einasto with 8 — 0.13) . We use 
Mollweide projection. To demonstrate the different morphologies, each template is normalized to 1 at the GC. Dark blue (dark 
grey) color refers to a flux suppressed by a factor of 10 -5 compared to the GC in each template. In eq.|8]we use the calculated 
flux values from the DRAGON package. 



members, a refers to the ratio of local DM density of 
the DD over the spherical halo given in eq. [3] The Iso^ 
template includes the contribution of DM extragalactic 
annihilations, the extragalactic background from other 
sources and the possible CR contamination. We also 
multiply the model map by the total FRONT+BACK- 
converted ULTRACLEAN photons exposure map and 
multiply by the mask. In Figure [3] for specific choices 
we show 4 different templates. On top left a specific 
model ("Back A") for the Galactic diffuse background 
is shown at the energies of interest. On top right we 
plot the " SphDM" template for an Einasto DM spherical 
halo profile with 8 = 0.13, and on bottom left we show 
the combined DM spherical and dark disk for a = 1 
(" SphDM" + "DarkDisk" + "MixedDM"). The DM sub- 
halos contribution ("SubDM") is given in the bottom 
right. We use HEALPix [3 with Nside = 128 which 
represents closer the angular resolution of Fermi LAT at 
these energies [do| . 



We also calculate the significance of a DM contribution 
from the diffuse by the test statistic, where 

^null 



TS = -2h 



(9) 



^bestfit 

£bestfit allows for the DM to contribute, while in C nu u 
we set the DM diffuse component to zero and refit the 
other two diffuse components. Our results are shown in 

Tables nnvi 

Ignoring first both the contribution of a dark disk and 
the DM subhalos, we find that the more cuspy DM pro- 
files for the main/spherical halo that lead to less DM 
contribution to the diffuse 7-ray spectrum away from the 
Galactic center, provide a larger positive fit to the An 
sky (see Table Yet the significance of that is not 
very large (TS = 9.1/9.2 and 11.9 in the best cases). An 
even more cored (fiat) Burkert DM profile for the main 
halo in the inner kpes does not give a better fit to the 
111 and 129 GeV lines distribution in agreement with the 
findings of and [l[ performed in subsections of the sky. 

Using the model of 8 = 0.13 for the Einasto DM den- 
sity profile that provides the best fit, we also test differ- 
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TABLE III: The values of relevant parameters for various assumptions on the Galactic diffuse background and the distribution 
of DM spherical halo ignoring the contribution of DD (a = 0) and subhalos. Second column gives the best fit annihilation 
cross-section assuming equal annihilation cross-sections to the two lines; {av)- 1 z = {<jv) 11 = av in units of xlO -28 cm 3 s _1 . 
Fi SO is the calculated isotropic flux at 111 ± 5 GeV and 129 ± 6 GeV in units of xl0~ 12 GeV -1 cm -2 s _1 sr . Columns 4-6 
refer to the 111 zb 5 GeV and 129 ± 6 GeV photons in the Background, DM and Isotropic diffuse components as predicted by 
the fitted values of N, A and B of eq. [S] Last column gives the TS for detection of DM signal. Values in parentheses refer to 
3a upper limits on the DM annihilation cross-section. 
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TABLE IV: The values of relevant parameters for various assumptions on the Galactic diffuse background, the distribution 
of DM spherical halo and the DD assuming the maximal DD contribution (a. = 1) and ignoring the contribution of subhalos. 
Second column is as in Table IHII Fi ao is in units of xl0~ 12 GeV -1 cm 2 s 1 sr 1 and av is in units of x 10 28 cm 3 s 1 . 
Columns 4-6 refer to the 111 ± 5 GeV and 129 ± 6 GeV photons in Back, DM and Iso components predicted by the fitted values 
of N, A and B. Last column gives the TS for detection of DM signal. 



ent diffuse background models ("Back A" -"Back D") to 
account for uncertainties in the interstellar medium gas 
distribution and interstellar radiation field distribution. 
We find that in all cases a DM component is preferred by 
the fit at ~2-3 a significance for the cuspier DM models 
(1-sided since in our fits we allow for the DM component 
to be even negative), and accounting for about 35-45 pho- 
tons (cross-sections of (av) = 1.3 — 2.7 x 10~ 28 cm 3 s _1 ) 
[661 ]. As with "Back A" the cuspier DM profiles are pre- 
ferred for the tested backgrounds. Our "Back A" is the 
same as the reference model described in [29( , which was 
shown to provide a good agreement to the 47r sky and in 
energies between 1 and 200 GeV and has also been cross- 
checked to local CR measurements. Model "Back B" as- 
sumes a different distribution for the molecular hydrogen 
gas component that is dominant at lower latitudes and 
toward the GC where many of the 111 and 129 GeV pho- 
ton excess have been claimed (see [1^] for more details). 
"Back C" and "Back D" Galactic diffuse models assume, 
respectively, an enhanced ISRF energy density distribu- 
tion toward the disk and a minimal metallicity gradient 
[52l |. The latter assumption affects both the morphology 
and the spectrum of the ISRF and as a result the 7-rays 
produced via Compton up-scattering of these photons by 
high energy CR electrons. These background models are 
discussed in further detail in (29|. They have not been 



constructed to study just the Galactic 7-ray background 
at ^110-130 GeV, but instead the general uncertainties 
in the Galactic diffuse 7-rays between 1-200 GeV, in the 
full sky and in subsections of it. 

Since in our fits we also allow for an isotropic com- 
ponent we can calculate the isotropic flux at these en- 
ergy ranges, taking into account also the Galactic diffuse 
background uncertainties. We find that ~ 210 — 230 pho- 
tons can be accounted by that component. This gives an 
isotropic flux of 5.6-6.1 xlO -12 GeV -1 cm~ 2 s _1 sr -1 
which is in agreement with the extragalactic isotropic 7- 
ray flux of [53| described by dN/dE cx £;- 2 - 41 ± 05 mea- 
sured between 200 MeV and 100 GeV. The ~ 210 - 230 
photons of the isotropic component can also be used to 
set approximate limits on the contribution from DM ha- 
los at extragalactic distances. 

We find that adding a dark disk component does not 
significantly change the fit to the data from the case of 
having only a spherical DM halo contribution. That is 
for any choice of diffuse background or spherical DM halo 
shown Table HVl Thus a significant DD contribution is 
not preferred by the current data; while yet it can not 
be excluded either. We give the likelihood fits for the 
maximal DD contribution. 
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TABLE V: The values of relevant parameters for various assumptions on the distribution of DM spherical halo, the DD and 
dim Galactic DM subhalos with different choices of the minimal subhalo mass. Units for av and Fi SO are as in Tables |IlI] and 
IIVI In parentheses we give the photon number associated to the subhalo diffuse component only. The subhalo mass spectral 
index is assumed to be -2. In all cases we use our reference "BackA" as Galactic diffuse background. 



Our fits give a relatively flat behavior in the parameter 
a with a slight preference for a — > 0. More specif- 
ically for a DD with zi/ 2 > 1.5 kpc (thick DD) with 
a — > 1 we get a lower test statistic by about ATS ~ 2 
compared to the cases of having only the spherical DM 
halos contributing, i.e. the data show a slight preference 
for having only a spherical DM halo at these energies. 
That is for any of our background assumptions. Yet thin 
DD (zi/2 — 0.5 kpc) give about the same TS as the only 
spherical DM halo contribution for any of our background 
models. Having tested also the 4 yr sample, that gave a 
slight preference for a = 1 with a zi/ 2 — 0.5 kpc DD 
we conclude that the current small number of photons 
at the two energy lines does not provide us with enough 
detail to discriminate a thin DD component from the 
disk-like distribution of the Galactic background 7-rays. 
With more 7-ray data we expect that we will be able to 
further disentangle the different 7-ray morphologies. 

Even for the case where we have a maximal DD contri- 
bution, i.e. when a = 1 most of the DM photons are due 
to the spherical and the mixed terms of eq. [SJ In those 
cases, the equivalent DM annihilation cross-section that 
accounts for the DM template normalization is maximal. 
That is the case since the total number of DM photons 

is about the same while the {^p 2 sph + p 2 DD + ^p S ph ■ Pdd^J 
term in eq. [S] is minimal. Since the DD term can not be 
excluded this is a way of allowing for higher annihilation 
cross-sections, by up to a factor of ~ 2. 

In Table IVl we also study the impact of the diffuse DM 
subhalo template " SubDM" of eq. [5] with and without 
a DD component. To bracket the uncertainties we used 
both the biased and the anti-biased Galactic subhalo dis- 
tributions (see Appendix). We also tested the subhalo 
mass extrapolations down to 10 +6 M Q and 1O _6 M . For 
the cases where the extrapolation is down to 10 +6 M Q 
the DM subhalo diffuse component is subdominant con- 
tributing only a few line photons. On the contrary, for ex- 
trapolations down to 10 _6 Mq that component can con- 



tribute up to - 1/3 (~ 1/2) of the Galactic DM line 
photons, with the remaining coming from the spherical 
main halo (spherical & DD components). With the cur- 
rent data the difference in the TS fit between the biased 
and the anti-biased distributions is very small (see Ta- 
ble [V}. Yet, even adding a strong subhalo term/template, 
our method can discriminate between different assump- 
tions for the cuspiness of the spherical halo " SphDM" . 
In agreement with all our previously discussed tests, a 
preference toward cuspier halos is found. In Table [V] we 
compare between 8 = 0.13 and S = 0.17 Einasto profiles 
that differ only in the inner few degrees from the GC. 

Having tested the DM case both with and without 
a dark disk and including/excluding the contribution of 
Galactic bound dim DM subhalos we have consistently 
found a thermally averaged cross-section (crv) = (av) 11 
= {<jv) 1 z for the two lines that is in the range of 1.5- 
4.5 xlO~ 28 cmV 1 . These values are a factor of 9-3 
smaller than the suggested values from analyses coming 
when concentrating only toward the GC Fitting the 
entire 4-7T sky can dilute the DM signal from the GC and 
thus suggest a smaller annihilation cross-section. Yet, we 
note that a strong annihilation annihilation cross-section 
to the lines should be seen at the diffuse spectrum at high 
latitudes when including the contribution from the dim 
subhalos with masses down to the free streaming scale. 
Our fits do not suggest such a case. Taking the subhalos 
to have masses down to 10 -6 M Q , we get from our fits 
shown in Table fVl ~ 15-20 photons (30 including a DD in 
the fit) from annihilations just in the dim subhalos. That 
number of photons from the diffuse subhalo component 
can be explained by a thermally averaged cross-section of 
1.5-2.5 xlO -28 cm 3 s _1 . Given that the velocity disper- 
sion in bound substructures is smaller than in the GC or 
locally that may be expected. Alternatively explained, 
thermally averaged cross-section of (av) = 1.0 x 10~ 27 
cm 3 s _1 down to the smallest subhalos would give ^100 
photons on the 47T sky just on the two lines and just from 
the dim Galactic subhalos. Of these ~ 25 photons would 
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be at | & | > 45° with the isotropic component predicting 
~ 70 photons in the two energy ranges (111 ± 5, 129 ± 6 
GeV). Thus the lines would have to be observed at high 
latitudes as well. If not found at these latitudes, then ei- 
ther the (av) is smaller in the subhalos, in general, or sub- 
halos are not formed down to masses of ~ 10" 6 -Mq. We 
note that in the above numbers we have not included the 
possible and more model dependent contribution form 
extragalactic DM annihilations. 



V. DISCUSSION AND CONCLUSIONS 

Recently, [l9[ has found indications for the 111 and 
the 129 GeV lines in the 2 yr Fermi unassociated point 
sources catalogue that would be indicative for DM anni- 
hilation in substructures. We compare the compatibility 
of the findings of [l9[ with the results from VLII cosmo- 
logical simulations and extrapolations of it. In that pro- 
cess we have assumed the same annihilation cross-section 
to the Galactic center lines signal [l-|y] given that the 
point sources lines signal and the Galactic center line(s) 
signal have been observed at the same 7-ray energies. We 
find that just considering VLII assumptions we do not 
get enough line photons from the brightest Galactic DM 
subhalos to claim agreement with the signal seen by [l!| . 
Yet completions of the VLII simulation results do give a 
number of line photons from the brightest subhalos that 
in the most optimistic cases is in good agreement with 
(but still below) the lines events number found by [19] (9 
events in their ULTRACLEAN sample) (see discussion in 
section irnj) . These most optimistic extrapolations to very 
small substructures (10 _6 M Q ) do not predict that many 
diffuse line photons to be conservatively constrained by 
the 7-ray flux at high latitudes. The same applies when 
we compare their predictions to the isotropic 7-ray flux 
which for the energies of the lines we have calculated to be 
centered at 5.6±0.3 xlO" 12 GeV" 1 cm" 2 s" 1 sr" 1 , with 
An sky fit values being in the range of 5.0 — 6.0 xlO" 12 
GeV" 1 cm" 2 s" 1 sr" 1 (accounting for systematic uncer- 
tainties in the evaluation of the isotropic component). 

Most of the 111 and the 129 GeV line photons on the 
7-ray sky are not from DM annihilations but rather from 
the Galactic diffuse and the isotropic diffuse background 
(shown in Fig. [2J . By doing a template fit to the Air sky, 
we test how robust the DM signal hypothesis is for dif- 
ferent physical assumptions on the Galactic diffuse 7-ray 
background flux and on the DM halo profile (see dis- 
cussion on section HV|) . We find a positive fit of having 
a DM spherical halo component. Our DM annihilation 
component hypothesis is preferred to the case without a 
DM component at a test statistic significance of up to 
12; with the values depending on the exact DM halo as- 
sumptions. More concentrated DM profiles give a larger 
significance to a DM signal. We also find that our results 
on the positive signal of DM annihilation are weakly de- 
pendent on Galactic diffuse background uncertainties re- 
lated to either the uncertainties in the distribution of the 



interstellar medium gas or to the energy density in the 
interstellar radiation field. 

Extending our set of tests on the DM distribution in 
the Galaxy, we include a dark disk component that could 
explain the non-isotropic distribution on the sky of the 
point sources of [19| . Our fits can not strongly favor 
or disfavor a significant disk- like DM component, even 
though there is a small preference toward thinner dark 
disks. We also study the contribution of the dimmest 
DM subhalos in the Milky Way to the diffuse gamma-ray 
sky at energies around 111 and 129 GeV and find a pref- 
erence toward an anti-biased distribution of the subhalos 
within the main spherical DM halo. Yet our analysis is 
somewhat constrained in its power by the small num- 
ber of photons on the sky at the energies of the 2 lines. 
As more statistics are being accumulated, a better un- 
derstanding of the morphology of the 7-ray sky will be 
achieved allowing for such a template analysis to further 
disentangle the background 7-ray sky from any possible 
DM component. 

Finally using the various combinations of backgrounds 
and DM distributions, we find the thermally averaged 
annihilation cross-section to be smaller than what has 
been originally suggested, with values ranging between 
1.5 - 4.5 x 10" 28 cm 3 s _1 . While the full sky fits are not 
optimal for a DM signal toward the GC, they include high 
latitudes which probe also the contribution from smaller 
substructures. A suppressed flux at high latitudes to 
the lines either indicates a smaller overall cross-section 
or a suppressed contribution from smaller substructures 
compared to the GC. 
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Appendix A: DM Signal from Biased and 
Anti-Biased Distributions of Subhalos in the Milky 

Way 

We fit the subhalos in the VLII simulation [|(| [54| with 
NFW profile [H: 

Psub(r) = Ps (jr) + • ( A1 ) 

and reading out the scale radius, r s , and scale density, 
p s , in terms of the peak of circular velocity V max and the 
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corresponding radius rymax (see, e.g. equation 9-10 in 
241), and their luminosities as: 



L = 4/3wripi l-(l + r t /r s y 



(A2) 



where r t is their tidal radius. By using this formula, 
combined with cq. ([5]), we find the results in the first 
row of Table U 

In deriving the biased distribution of subhalos, first, 
we calculate the gravitational potential from the host 
halo, cutting the profile at its virial radius, and distin- 
guish bound subhalos from the unbound ones. There are 
13510 bound subhalos, 9372 of them are within virial ra- 
dius. We also calculate each subhalo's pericenter, r p , 
to obtain the pericenter probability distribution func- 
tion (PDF) as a function of Galactocentric radius, r g , e.g 
dP/dr p (r p , r g ). For unbound subhalos, we identify their 
tidal radius as their virial radius, and their tidal mass 
as their virial mass. According to virial theorem, the 
average density inside virial radius should be the same 
for all halos. We find that this is indeed the case in 
VLII simulation: the average density inside virial ra- 
dius of unbound subhalos is Aop cr it, where Ao is the 
overdensity relative to critical matter density for spher- 
ical collapse for z = 0, and p cr it is the critical mat- 
ter density. For bound subhalos, tidal radius is usually 
much smaller than virial radius because of tidal strip- 
ping. We define the tidal concentration of a subhalo, 
Ct, as its average density within tidal radius divided by 
critical matter density; for unbound subhalos, Ct = Aq. 
To determine the virial concentration PDF for unbound 
subhalos, we use the Bullock model [55[ as parametrized 
in eq. (2) of [56] which relates the median virial con- 
centration of a subhalo, c V i r — r V i r /r s , with its virial 
mass, m v i r . We find that this parametrization can fit 
unbound subhalos's virial concentrations fairly well. We 
take the mass PDF of unbound subhalos to be a power 
law, dP/dm t oc m^~ a O(m t — m cut ), and vary the spectral 
index a = [1.9,2]. In fitting the mass and concentration 
PDF, we use the unbound subhalos because they are not 
tidally stripped. Hence, in this model, the Galactocentric 
radial dependence of mass and concentration will appear 
later, after they are tidally stripped by their host. To 
find subhalo's minimum concentration at some Galacto- 
centric radius r g , we apply the Roche criteria: for a sub- 
halo in circular orbit, subhalo's self-gravity at r t should 
be equal to the differential gravity pull of the host halo 
computed at r g . As the subhalo's orbit is not exactly 
circular, the tidal force is strongest at its pericenter, so 
the concentration is determined by the tidal forces at its 
pericenter [571 ]: 



Ctr{r p ) 



2ph{< r p ) 3ph{r P ) 



Pcrit 



Pcrit 



(A3) 



where ph is VLII host density profile. We then refine the 
previous estimate by taking into account the pericenter 
radius distribution for subhalos, dP/dr p . Subhalos ini- 
tially with ct = Ao > ctr {r p ) are left intact, whereas 



subhalos initially with ct < Ct r (r p ) are tidally stripped 
until c tr ( r p) is reached, e.g.: 



c[ = max [c tr (r p ) , A ] 



(A4) 



It is a fair approximation that scale p s and r s do not 
change in this process 58]. By following these steps; 
a subhalo's mass, concentration, and luminosity after 
tidal stripping are completely determined by its pericen- 
ter, mass, and concentration before tidal stripping. We 
then calculate the average mass, concentration, and lu- 
minosity after tidal stripping, and compare them against 
simulation. We find that luminosity and mass can be 
fitted very well, with best fit parameters {m cut ,a) — 
(10 6 - 5 , 1.9). (The m cut here is the minimum subhalo 
mass before tidal stripping.) Our model also fits the 
tidal concentration very well near host center, although 
slightly deviates near its virial radius. This might indi- 
cate that our treatment of dynamical effects is oversimpli- 
fied. Especially, the minimum concentration calculated 
by Rochi criteria is only achieved after several pericentric 
passages (5m60| . whereas in our model, we assume that 
all of them are already above minimum concentration. 
However, the procedure is validated by very good fits to 
the luminosities. Besides, it is conceptually simple and 
can be easily generalized to other host and subhalo mass 
profiles, and redshifts. 

Regarding the subhalos number density profile, we con- 
sider both biased and anti-biased case. As mentioned 
in the section IIII1 the strong tidal force near GC could 
make the subhalos spatial distribution to be anti-biased 
with respect to the host density profile. Indeed, for VLII 
itself, we find that the deviation from biased starts to 
happen around ~ 30 kpc. However, we find that the sub- 
halo spatial distribution for VLII is more biased with 
respect to VLI [13, [H[ , which has lower resolution. This 
might indicate that the anti-biased is the result of nu- 
merical effect (The strong tidal force could strip subhalos 
until they are below or near resolution limit, hence un- 
detectable as subhalos.) Also, VLII subhalos which are 
selected by mass show more anti-biased tendency than 
the ones which are selected by maximum circular veloc- 
ity, confirming [23|,[4Cj. Keeping these in mind, we take 
biased distribution for subhalos number density: 



n 



* j sub 



= 5.84 kpc 



m cu t 
,10 6 - 5 M ( 



-a+l 



Ph (r) 
M Q pc- 3 



(A5) 



With this normalization, the total mass and number 
of subhalos with tidal mass bigger than 10 6 5 between 
~ 30 kpc and r" r in VLII simulation is within 5% from 
the value calculated with our model using best fit param- 
eters. 

Aquarius simulation Aq-A-1 [4lj is another simulation 
which has parameters similar to VLII. Specifically, it has 
nominal mass resolution, m p = 1.71 x 10 3 M©, host halo 
mass, M 50 = 2.5 23 x 10 12 M Q , and host halo radius, r 50 = 
433 kpc. However, their number of subhalos with mass 
bigger than ~ 10 7 M Q is approximately twice VLB's. We 
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discuss the modification to our results in Table Q] and 
Table [TT] if we use Aquarius normalization. 

From this procedure, we can find the number den- 
sity of subhalos per unit luminosity after tidal stripping, 
dn su i,/dL. The number of photons that we receive from 
a single subhalo with luminosity L and line of sight (los) 
distance A from us, for channel ch = 77 or 7Z is given 
by (0. By folding ([5]) with dn su h/dL, we can obtain the 
values in Table Q] and Table [TTJ On the other hand, the 
contribution to the host's (p 2 ) from subhalos is given by: 

(pLb) = n sub x (L) , (A6) 

where ( L) is the average luminosity of the entire subhalos 
population after tidal stripping. 

The anti-biased distribution is taken from appendix 
A of [HI- The normalization is such that the total mass 
in subhalos with masses between \Q~ 5 Mh and 10 _2 M^ is 
3.4% of Mh. (In the original paper, the authors normalize 
to 10% of Mj,. However, if we only include bound subhalos 



within virial radius, the mass fraction is only 3.4%.) In 
this model, the authors encapsulates the tidal force from 
the host halo by adding radial dependence to the m vir — 
c V i r relation, so that subhalo closer to GC has higher 
concentration. 

As a further refinement, we also add the Galactic disk 
to the VLII halo. For the Galactic disk model, we fol- 
low [35[, for NFW parameter and our distance to GC, 
Rq = 8.5 kpc. We spherically average the Galactic disk. 
The relative difference in density between the averaged 
version and the original version is substantial only in the 
disk plane. There is no noticeable difference to the peri- 
center PDF, and the only modification to biased distri- 
bution is that, we replace ph — > Ph + Pdisk in (|A3|1 . For 
anti-biased one, the effect of host tidal forces have been 
taken into account by adding radial dependence to the 
c v i r — m v i r relation. Therefore, in eq. (|A3|) . we replace 
Ph Pdisk- We also use the same pericenter distribution 
as in the biased one. 



C. Weniger (2012), 1204.2797. 

M. Su and D. P. Finkbeiner (2012), 1206.1616. 

E. Tempel, A. Hektor, and M. Raidal (2012), 1205.1045. 

A. Rajaraman, T. M. Tait, and D. Whiteson (2012), 

1205.4723. 

W. Buchmuller and M. Garny (2012), 1206.7056. 

A. Boyarsky, D. Malyshev, and O. Ruchayskiy (2012), 

1205.4700. 

M. Ackermann et al. (LAT Collaboration) (2012), 
1205.2739. 

M. Kuhlen, J. Guedes, A. Pillepich, P. Madau, and 
L. Mayer (2012), 1208.4844. 
K. Rao and D. Whiteson (2012), 1210.4934. 
E. Charles, A. Albert, E. Bloom, S. A., Z. Yang, and 
G. Zaharijas, J^th fermi symposium, monterey ca (2012). 
A. W. B. Albert, H. R., and E. Bloom, Jfih fermi sympo- 
sium, monterey ca (2012). 

D. P. Finkbeiner, M. Su, 
1209.4562. 

D. Whiteson, JCAP 1211, 008 (2012), 
D. Whiteson (2013), 1302.0427. 
M. R. Buckley and D. Hooper (2012), 1205.6811. 
T. Cohen, M. Lisanti, T. R. Slatyer, and J. G. Wacker 
(2012), 1207.0800. 

I. Cholis, M. Tavakoli, and P. Ullio, Phys.Rev. D86, 
083525 (2012), 1207.1468. 

A. Hektor, M. Raidal, and E. Tempel (2012), 1207.4466. 
M. Su and D. P. Finkbeiner (2012), 1207.7060. 
A. Hektor, M. Raidal, and E. Tempel (2012), 1208.1996. 
D. Hooper and T. Linden, Phys.Rev. D86, 083532 
(2012), 1208.0828. 

J. Diemand, M. Kuhlen, and P. Madau, Astrophys.J. 
657, 262 (2007), astro-ph/0611370. 

J. Diemand, M. Kuhlen, and P. Madau, Astrophys.J. 
667, 859 (2007), astro-ph/0703337. 

M. Kuhlen, J. Diemand, and P. Madau, Astrophys.J. 
686, 262 (2008), 0805.4416. 

C. Evoli, D. Gaggero, D. Grasso, and L. Maccione, JCAP 



and C. Weniger (2012), 



1208.3677. 



[26 
[27 

[28 
[29 



[30 
[31 
[32 
[33 
[34 
[35 
[36 
[37 
[38 

[39 

[40 
[41 

[42 

[43 
[44 



0810, 018 (2008), 0807.4730. 
http: / /www. desy.de / maccione/DRAGON/. 
I. Cholis, M. Tavakoli, C. Evoli, L. Maccione, and P. Ul- 
lio, JCAP 1205, 004 (2012), 1106.5073. 
M. Tavakoli (2012), 1207.6150. 

M. Tavakoli, I. Cholis, C. Evoli, and P. Ullio (2013), Dif- 
fuse Galactic Gamma Rays at intermediate and high lat- 
itudes. II. Constraints on the DM properties. 
I. A. Grenier, J.-M. Casandjian, and R. Terrier, Science 
307, 1292 (2005). 

T. F.-L. Collaboration (The Fermi-LAT Collaboration), 

Astrophys.J. 750, 3 (2012), 1202.4039. 

T. Bringmann, X. Huang, A. Ibarra, S. Vogl, and 

C. Weniger (2012), 1203.1312. 

D. Merritt, J. F. Navarro, A. Ludlow, and A. Jenkins, 
Astrophys. J. Lett. 624, L85 (2005), astro-ph/0502515. 
J. I. Read, G. Lake, O. Agertz, and V. P. Debattista 
(2008), 0803.2714. 
R. Catena and P. 
0907.0018. 

P. Salucci, F. Nesti, 
tron. Astrophys. 523, 
P. M. W. Kalberla, 

(2007) , 0704.3925. 
A. M. Green, S. 
Mon.Not.Roy.Astron.Soc. 
ph/0309621. 

S. Profumo, K. Sigurdson, and M. Kamionkowski, 
Phys.Rev.Lett. 97, 031301 (2006), astro-ph/0603373. 
J. Diemand et al., Nature 454, 735 (2008), 0805.1244. 
V. Springel, J. Wang, M. Vogelsberger, A. Ludlow, 
A. Jenkins, et al., Mon.Not.Roy.Astron.Soc. 391, 1685 

(2008) , 0809.0898. 

P. L. Nolan et al. (Fermi-LAT Collaboration), Astro- 
phys.J. Suppl. 199, 31 (2012), 1108.1435. 
A. Drlica- Wagner, Private communication (2012). 
K. Gorski, E. Hivon, A. Banday, B. Wandelt, F. Hansen, 
et al., Astrophys.J. 622, 759 (2005), astro-ph/0409513. 



Ullio, JCAP 1008, 004 (2010), 

G. Gentile, and C. Martins, As- 
A83 (2010), 1003.3101. 
L. Dedes, J. Kerp, and U. Haud 



Hofmann, 
353, 



and D. J. Schwarz, 
L23 (2004), astro- 



13 



[45] N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer, and 

N. Weiner, Phys.Rev. D79, 015014 (2009), 0810.0713. 
[46] M. Lattanzi and J. I. Silk, Phys.Rev. D79, 083523 

(2009), 0812.0360. 
[47] J. Hisano, S. Matsumoto, and M. M. Nojiri, Phys. Rev. 

Lett. 92, 031303 (2004), hep-ph/0307216. 
[48] J. Hisano, S. Matsumoto, M. M. Nojiri, and O. Saito, 

Phys. Rev. DTI, 063528 (2005), hep-ph/0412403. 
[49] I. Cholis and L. Goodenough, JCAP 1009, 010 (2010), 

1006.2089. 

[50] T. R. Slatyer, N. Toro, and N. Weiner, Phys.Rev. D86, 
083534 (2012), 1107.3546. 

[51] G. Dobler, D. P. Finkbeiner, I. Cholis, T. R. Slatyer, and 
N. Weiner, Astrophys.J. 717, 825 (2010), 0910.4583. 

[52] T. A. Porter, I. V. Moskalenko, A. W. Strong, E. Or- 
lando, and L. Bouchet, Astrophys. J. 682, 400 (2008), 
0804.1774. 

[53] A. Abdo et al. (The Fermi-LAT collaboration), 

Phys.Rev.Lett. 104, 101101 (2010), 1002.3603. 
[54] http:/ /www. ucolick.org/~diemand/vl/data.html. 
[55] J. S. Bullock, T. S. Kolatt, Y. Sigad, R. S. Somerville, 

A. V. Kravtsov, et al., Mon.Not.Roy.Astron.Soc. 321, 

559 (2001), astro-ph/9908159. 
[56] H. Zechlin, M. Fernandes, D. Elsaesser, and D. Horns, 

Astron.Astrophys. 538, A93 (2012), 1111.3514. 
[57] G. Tormen, A. Diaferio, and D. Syer (1997), astro- 

ph/9712222. 

[58] S. Kazantzidis, L. Mayer, C. Mastropietro, J. Diemand, 
J. Stadel, et al., Astrophys.J. 608, 663 (2004), astro- 
ph/0312194. 

[59] E. Hayashi, J. F. Navarro, J. E. Taylor, J. Stadel, 
and T. R. Quinn, Astrophys.J. 584, 541 (2003), astro- 
ph/0203004. 

[60] J. Gan, X. Kang, F. C. den Bosch, and J. Hou (2010), 
1007.0023. 



[61] There have beenl6 unassociated point sources in the 2yr 
Fermi catalogue detected at the level of ~4cr or more. In 
oder to reach that significance a spectrum at low energies 
is also needed. DM does not give a strong spectrum signal 
at low energies so the actual spectrum of these sources 
(if they are DM subhalos) may be contaminated by other 
close- by point sources or diffuse 7-rays. That may also 
result in a non-trivial selection effect. 

[62] Our 4.4 yr sample refers to the data from 04 August 
2008-03 January 2013. 

[63] This value is obtained assuming isotropic distribution. 
The actual distribution of the 111 and 129 GeV photons 
is not isotropic but rather peaks on the Galactic disk. 
Our results do not differ that much though since the dis- 
tribution of the 575 unassociated point sources also peaks 
on the Galactic disk. 

[64] fl9| suggest a p-value of p = 0.00069 for having 13 
SOURCE events given a background described by a 
power-law. 

[65] The angular resolution of the Fermi LAT instrument is 
different at the energies the we study by about a factor 
of 2 between the front and the back-converted converted 
photons. With Nside=128 our pixels include entirely the 
front-converted photons are slightly smaller than the 68% 
containment angle for back-converted events. Choosing 
Nside=256 would have resulted in photon events being 
spread in multiple pixels thus making them correlated; 
while Nside< 64, would be underestimating the ability 
of the instrument to resolve 7-ray structures at these 
energies. We have checked that our best fit results are 
consistent within 2% for Nside of 128, 64, 32, 16 and 8. 

[66] Since we use the entire 4tt sky and have selected in ad- 
vance the energy range of the 7-ray events based on the 
works of there is no look elsewhere effect accounted 
for. 



